Accepted for publication in the Astrophysical Journal 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



ISOTROPIC MID-INFRARED EMISSION FROM THE CENTRAL 100 PC OF ACTIVE GALAXIES 
N. A. Levenson 1,2 , J. T. Radomski 2 , C. Packham 3 , R. E. Mason 4 , J. J. Schaefer 3 , and C. M. Telesco 3 



On 
O 
O 
(N 

bJQ' 
< 



O 

u 

6 

c3 



> 

m 
o 
q 

o 

OS 

o 



- 1—^ 

x 



ABSTRACT 

Dust reprocesses the intrinsic radiation of active galactic nuclei (AGNs) to emerge at longer wave- 
lengths. The observed mid-infrared (MIR) luminosity depends fundamentally on the luminosity of the 
central engine, but in detail it also depends on the geometric distribution of the surrounding dust. To 
quantify this relationship, we observe nearby normal AGNs in the MIR to achieve spatial resolution 
better than 100 pc, and we use absorption-corrected X-ray luminosity as a proxy for the intrinsic AGN 
emission. We find no significant difference between optically classified Seyfert 1 and 2 galaxies. Spec- 
troscopic differences, both at optical and IR wavelengths, indicate that the immediate surroundings 
of AGNs is not spherically symmetric, as in standard unified AGN models. A quantitative analysis 
of clumpy torus radiative transfer models shows that a clumpy local environment can account for 
this dependence on viewing geometry while producing MIR continuum emission that remains nearly 
isotropic, as we observe, although the material is not optically thin at these wavelengths. We find 
some luminosity dependence on the X-ray/MIR correlation in the smallest scale measurements, which 
may indicate enhanced dust emission associated with star formation, even on these sub-100 pc scales. 
Subject headings: galaxies: active — galaxies: Seyfert — infrared: galaxies 



1. INTRODUCTION 

While active galactic nuclei (AGNs) present a variety 
of observational characteristics, unified models suggest 
that the central engines of all these objects are fundamen- 
tally the same. Accretion onto a supermassive black hole 
is the basic energy source, and an optically and geomet- 
rically thick surrounding "torus" introduces effects that 
depend on viewing geometry. In particular, spectrally 
broad emission lines are observed in type 1 AGNs, which 
afford a direct view of the fast-moving material close to 
the central engine, while the torus b locks the view of t his 
broad line region in type 2 objects (|Antonucci 1119931 ). 

The dust in the obscuring torus also reprocesses the in- 
trinsic hard continuum radiation to longer wavelengths, 
and the bulk of it emerges in the mid-infrared (MIR) 
regime (5-30/im). The observed MIR emission of AGNs 
thus is sensitive to the intrinsic bolometric luminosity of 
the AGN. High energy X-ray emission reveals the bolo- 
metric luminosity of the AGNs, and few other energy 
sources contribute substantially to these observed fluxes. 
Provided that measured 2-10 keV luminosities are cor- 
rected for intrinsic absorption, which can be significant, 
this X-ray luminosity is an effective proxy for the AGN 
bolo metric luminosity, representing approximately 5% of 
L bol (lElvis et al.l ll99l. 

In detail, the spectral energy distribution of the 
reprocessed emission also depends on the geometric 
distribution of the dust. If the torus material is 
smoothly distributed, for example, the observed IR flux 
along unobscured lines of sight to the hot, optically- 
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thin surface of the inner torus is much greater than 
that measured along lines of sight through the opti- 
cally thick torus dPier fc Kro lik 1992; lGranato fc Danesd 
119941: lEfstathiou fc Rowan-Robinsonl(l995h . In contrast, 
if the dusty material is separated into clumps that do not 
fill the torus volume, t he luminosity dependence on view- 
ing angle is reduced (iHonig et al.l 120061 : iNenkova et al.l 
I2008bt ISchartmann et alj 120081 ). MIR photometry of 
type 1 and 2 AGNs can therefore sensitively discrim- 
inate between smooth and clumpy dust distributions. 
Moreover, recent observations at high spatial resolution 
(iJaffe et al.1 l2004t iPackham et al.1 120051: iTrist ram et all 
l2007tlRadomski et al.l2008h suggest that the spatial scale 
of the dusty torus is small (< 10 pc). Thus, measure- 
ments of the torus emission must be made over limited 
physical regions, and comparison models must accommo- 
date the small torus extent. 

Previous work generally indicates that the intrinsic 
and reprocessed AGN emission are strongly correlated. 
However, s ome of thes e stud ies su ffer from poor angular 
resolu tion. iLutz et ail (|2004f ) and lRamos Almeida et al.l 
(|2007ft . for example, use MIR data from ISO. The spa- 
tial resolution of these observations is typically 2 kpc, 
and dust emission associated with star formation can 
contribute greatly on these scales. Observations with 
8m class telescopes on the ground offer an order of 
magnitude greater angular resolution, but initial in- 
vestigations provi ded only small samples (|Krabbe et al.l 
| 200lf iHorst et al.l I2006T). Extending the samples (e.g., 
iHorst et all 120081 : iGandhi et al.ll2009r i introduced mixed 
classes of AGNs (including radio galaxies and LINERs, 
for example) which may have distinct properties because 
of their radio loudness and hardness of the ionizing con- 
tinuum. Despite these complications, no earlier work 
shows significant differences between type 1 and 2 AGNs 
in the MIR. The spectral energy distributions of type 1 
and 2 AGNs are characteristically different through the 
near-IR, with the former emitting much more strongly. 
The only suggestion of a difference between types in the 
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MIR emerges at relatively sh ort wavelengths (A = 6.7/im; 
iRamos Almeida etP al. 2007) and is not maintained in the 
longer-wavelength data of the same study (C. Ramos 
Almeida, private communication 2009). 

In this paper we present MIR imaging of a sample 
of 17 Seyfert galaxies to compare with X-ray measure- 
ments of the AGN bolometric luminosity. We restrict 
this study to normal Seyfert galaxies, excluding radio- 
loud sources, in which non-thermal synchrotron emission 
can contribute significantly to the MIR flux. We require 
the sample galaxies to be closer than 50 Mpc, to achieve 
spatial resolution around 100 pc. In contrast with other 
work, we do not consider LI NER galaxies, which may 
have unusual M IR properties (|Elitzur fc Shlosm"anll2006t 
iPerhnan et al.ll2~001l ). Moreover, we take particular care 
to correct the observed X-ray fluxes to recover the in- 
trinsic AGN emission. The X-rays advantageously offer 
observations of emission that originates close to the black 
hole, avoiding the uncertainties of extinction corrections 
and variations among the more distant narrow line re- 
gions of optical [O III] measurements, for example. In 
addition, we directly compare our MIR results with ra- 
diative transfer calculations of dust emission. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. MIR Measurements 



Observations were made using T-ReCS dTelesco et al 



1998) on the 8.1m Gemini South, Michelle (Glass e et al 



19971) on Gemini North, and OSCIR on the Blanco 4m 
at CTIO, Gemini North, and Gemini South. Table [Q 
shows the log of observations. OSCIR used a Boeing 
128 x 128 pixel Si:As BIB detector, providing a plate 
scale of 0.183" per pixel on the 4m at CTIO. The spa- 
tial resolution achieved was around 1". Both T-ReCS 
and Michelle use a Raytheon 320 x 240 pixel Si:As IBC 
array, providing a plate scale of 0.089 and 0.1005" per 
pixel respectively, and achieved a spatial resolution of 
between 0.3 to 0.5". In the case of T-ReCS, the detector 
was used in corre lated quadruple sampling (CQS) mode 
dSako et al.ll200l . In all cases, images were obtained 
using the standard chop-nod technique to remove time- 
variable sky background, telescope thermal emission and 
so-called "1//" detector noise. For the T-ReCS obser- 
vations, the chop throw was 15" and the telescope was 
nodded every 30 s, whereas Michelle was nodded every 
45 s. The OSCIR observations had a chop throw of 30" 
on the 4m, but 15" on Gemini, and all telescopes were 
nodded every 30 s. The chop throw was typically at 0° 
(N-S), but was angled in the case of extended objects. 
Observations were made using the N (A c = 10.36/im, 
AA = 5.27^m), N' (A c = 11.2/im, AA = 2.4/xm), Si2 
(A c = 8.74/im, AA = 0.78/im), or Si5 (A c = 11.6/xm, 
AA = 1.1/zm) filters, where A c is the central wavelength, 
and the filter width AA indicates 50% cut-on and cut-off 
values. 

Data were reduced using in-house developed IDL or 
equivalent Gemini IRAF routines. The difference for 
each chopped pair for each given nod-set was calculated, 
and the nod-sets were then differenced and combined to 
create a single image. Chopped pairs obviously com- 
promised by cirrus, high electronic noise, or other prob- 
lems were excluded. OSCIR, T-ReCS, and Michelle were 
mounted on the Cassegrain port of the telescopes so that 
north was up and east was left as projected onto the 



detector. In post-processing, images of the PSF stars 
were de-rotated to match the telescope pupil PA when 
the galaxy observations were observed. Image rotation is 
necessary because the projected angles of the telescope 
pupil (particularly the secondary mirror supports) rotate 
during observations (or pointings), having a significant 
effect on the low-level profile. 

Both flux and point source function (PSF) standard 
observations were made for flux and image quality cal- 
ibrations through the same filters used for each galaxy 
observation. Flux standards were observed using stan- 
dard sources suggested from the Gemini web pages. We 
find that, consistent with other mid-IR observations, the 
uncertainty in the flux determination is around 10%. 
PSF observations were made using the same filters as 
the galaxy observations, immediately prior to or after 
the galaxy observations and used an identical setup to 
accurately sample the image quality. Short PSF or flux 
standard observations are comparable to longer source 
observations as the closed-loop active optics of Gemini 
and the mechanical stability offered by the 4m at the 
CTIO provides a similar PSF when taken at a similar 
telescope pointing and time. Observations of PSF and 
flux standards through the night showed a stable and 
consistent PSF. 

We made three different photometric measurements 
of each nucleus. First, assuming that the peak emis- 
sion represents the central source, we scaled the PSF 
to this peak. The scaled flux integrated over a region 
that fully encompasses the Airy disk is the "diffraction- 
limited" measurement, where diameters of 2" are typi- 
cal of the observations from the 8.1m Gemini telescopes 
and 4" is appropriate for the observations from the 4m 
Blanco telescope. Second, we used a fixed physical scale 
defined by the resolution of the most distant galaxy, of 
100 pc. Third, we estimated the flux of the unresolved 
nuclear component, the "PSF-fitting" measurement. In 
this case, we subtracted a scaled PSF from the galaxy 
image. The PSF scaled to the peak of the galaxy emis- 
sion (the "diffraction-limited" measurement) represents 
the maximum contribution of the unresolved source. The 
residual of the total emission minus the scaled PSF rep- 
resents the host galaxy contribution. [A detailed study 
of the extended near nuclear structures of the galaxies 
in our sample will be the subject of a forthcoming pa- 
per (J. T. Radomski et al., in preparation).] Subtrac- 
tion of a PSF scaled to match the peak of the galaxy 
emission exactly results in an unrealistic minimum in 
the residual host galaxy emission. Instead, we reduce 
the scale of the PSF to produce a "flat" profile in the 
residual to obtain the unresolved fluxes reported in Ta- 
ble [2 (See[R amos Almeida et al.ll2009l for examples of 
this technique.) We note that a flat nuclear profile may 
not account for any residual stellar cusp at the nucleus, 
but this effect will be small. For reproducibility and sim- 
plicity, we apply this approach to all of our sources. We 
estimate this method introduces a further uncertainty of 
10-15% in the final flux measurements, which we add 
in quadrature to the flux calibration errors. In the case 
of the fixed spatial scale, a simple aperture was placed 
over the galaxy and the flux extracted. In all cases, a 
sky annulus was used to subtract low-level, residual sky 
emission not cancelled by the chopping and nodding. 

While the observations were made through different 
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filters in the 10/zm window, we assume a flat spectral en- 
ergy distribution through the limited wavelength range 
they cover. Thus, we make no corrections to the mea- 
sured fluxes in the comparisons below. We expect this 
simplified approach minimizes the introduction of addi- 
tional uncertainties, given that we do not know in detail 
the high spatial resolution spectral energy distribution 
of each AGN. We have few of our own observations of 
Seyfert 1 galaxies, so we considered supplementing these 
with results from the literature. However, in order to 
compare the different MIR measurement techniques, we 
quote the statistical results from our sample alone. 

2.2. X-ray Measurements 

We take care that the X-ray measurements also isolate 
the AGN, favoring high angular resolution data. We pre- 
fer observations close in time to the MIR observations, to 
avoid possible complications of intrinsic AGN variability. 
All the X-ray fluxes are corrected for line of sight absorp- 
tion (due to the host galaxy, or Milky Way, or both) to 
determine intrinsic luminosity in the 2-10 keV band, Lx- 
The underlying AGN continuum of Compton thick galax- 
ies, those obscured by Nji > 10 24 cm' 2 , does not emerge 
directly below 10 keV. In these cases, we prefer higher- 
energy observations to determine the AGN luminosity. 
These are generally uncontaminated by stellar processes, 
despite the poor angular resolution. 

Table [5] lists the sample and basic data. The tabulated 
uncertainties reflect the range of observed luminosities, 
and genuine variability is significant in Seyfert Is. The 
span of measured Seyfert 2 luminosities is usually small 
for each AGN, and we adopt a nominal 20% uncertainty 
to account for intrinsic variability if none is observed 
directly. When detection of a Compton thick AGN at 
E > 10 keV is absent, the tabulated range indicates the 
uncertainty in the correction to recover the intrinsic lumi- 
nosity. We use measurements from the literature where 
they are available. We present our own X-ray analysis of 
three galaxies below, and we give special consideration 
to the exceptional case of NGC 1068. 

NGC 1068 is extremely optically thick, and the di- 
rect AGN emissi on is not detected even above 10 keV 
(|Matt et al.lll997l ). In this case, the range of Lx indicates 
the uncertainty in det ermining the int r insic luminosity. 
We use the results of iLevenson et al.l ([2006) , adopting 
the estimate of Lx from the Fe Kq line equivalent width 
(EW) and luminosity, and considering the fitted reflec- 
tion model a l ower limit, for D = 14.4 Mpc. We use 
the estimate of [Iwasawa et al.l (|l997), who consider the 
properties of the ionized scattering region, as an upper 
limit on Lx- 

NGC 1566 is included in the XMM slew survey and 
detec ted in both hard a nd soft bands in one observa- 
tion (|Saxton et al.ll2008h . Assuming T = 1.7 (to ap- 
ply the standard conversion from the survey count rate 
to flux) we have absorption-corrected 2-10 keV flux 
F x = 5.8(±2) x 10" 12 erg ernes' 1 . 

We fit the ASCA spectrum of NGC 3081, using source 
and background screened events files from the Tartarus 
database. We simultaneously fit the 3-10 keV data from 
all four detectors, allowing for a constant offset due to 
calibration uncertainty. We ignore the softest energies 
because they do not directly reveal the intrinsic AGN 
output. Instead, they show significant line emission, 



suggesting a strong thermal or photoionized contribu- 
tion. We fit the AGN as an absorbed power law, fix- 
ing photon index T = 1.9, which is typical of Seyfert 
galaxies. We find N H = 6.3(±0.4) x 10 23 cm" 2 , with 
absorption-corrected 2-10 keV flux F x = 3.8(±0.2) x 
10 -11 ergcm~ 2 s _1 averaged over all detectors. 

We fit the Chandra spectrum of the nucleus of 
NGC 5728. The da ta were reduced as described in 
ILevenson et al.l (|2006f ). and we used a 2" nuclear aper- 
ture. The direct continuum emission is strong. Fit- 
ting the data at E > 3keV and fixing T = 1.9, we 
find N H = 7.1(±2) x 10 23 cm" 2 , and F x = 9.3(±2) x 
10 -12 ergcm~ 2 s _1 . The neutral Fe Ka line is also 
strong, having EW = 790 eV. 

3. RESULTS AND DISCUSSION 

3.1. MIR and Bolometric Luminosity Correlations 

Consistent with earlier similar work (|Horst et al.ll2006l 
l2008UGandhi et al . 1 120091 ) . we find overall that the X-ray 
and MIR luminosities are well-correlated. Because the 
X-ray measurements represent the intrinsic AGN lumi- 
nosity, with minimal contamination from other sources 
in the host galaxies, we conclude that the MIR emis- 
sion is strongly correlated with the AGN bolometric lu- 
minosity. This result holds for all three methods we 
employ to measure the MIR emission, and considering 
both intrinsic variability of the X-ray emission and using 
only errors on individual X-ray measurements. One sim- 
ple way to quantify the strong correlation is to measure 
the unweighted mean (/z) and standard deviation (a r ) of 
log Lm i r/Lx- Using the diffraction-limited MIR lumi- 
nosity we find /i = 0.12 and a r = 0.6; in the PSF-fitted 
case we find n = 0.10 and a r — 0.7, and with the fixed 
physical aperture we find fi — 0.37 and oy = 0.5. These 
values are comparable to the results of I Gandhi et al.l 
(2009), although we measure greater dispersion in the 
present smaller sample. 

Despite the intrinsic high resolution of the diffraction- 
limited measurements, the differences with the PSF- 
fitted fluxes indicates that even on large telescopes, 
diffraction-limited MIR observations do not truly isolate 
the immediate dusty surroundings of AGNs, which ar e 
confined to smaller (< 10 pc) scales (jMason et al.| [2006'). 
Given that we cannot measure the torus by itself, the 
problem of Lf IIR is that the physical scale it encom- 
passes varies among the different galaxies. Thus, while it 
offers a consistent technique from the observational view- 
point, the physical correspondence changes. The use of 
a fixed physical aperture directly compensates for this 
disadvantage, while reducing the effective resolution (for 
all but one galaxy). The PSF fitting offers the best de- 
termination of the flux of the unresolved MIR-emitting 
torus, despite the variation of physical scale from galaxy 
to galaxy. Thus, this result best describes the relation- 
ship between intrinsic AGN luminosity and the repro- 
cessed emission of the torus. 

The MIR luminosity of NGC 4945 changes greatly with 
the three different measurement techniques. This galaxy 
contains a strong starburst. Other sources of MIR emis- 
sion in addition to the AGN are therefore significant on 
small scales, and these dominate even Lf IIR , unlike in 
other galaxies. Also, because NGC 4945 is nearby, the 
fixed 100-pc aperture measurement covers a much larger 
area (that contains strong emission sources) than the 
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diffraction-limited (36 pc) region does. 

The differences among the various MIR measurements 
are larger for the Seyfert 2 galaxies than for the Seyfert 
Is. Among the latter, the PSF-fitting scales are all 
90 or 100%, producing nearly the same result as the 
diffraction-limited measurements. The fixed physical 
aperture measurements remain similar for the Seyfert Is. 
In contrast, PSF fitting of some Seyfert 2s implies that 
only a small fraction of even the unresolved emission is 
the torus. In these galaxies, other contributions remain 
significant on the 100 pc scale, as well. 

Table [3] lists the Spearman rank correlation coefficient 
and the corresponding probability, considering Seyfert 
1, Seyfert 2, and total samples, and the three differ- 
ent MIR measurements separately. The correlations of 
the Seyfert 2 and full samples are significant; the less 
significant (though strong) correlation of the Seyfert 1 
subsample alone is a consequence of the small number 
of objects. Figures [T] through [3] show the data and fit- 
ted functions, where X-ray variability contributes to the 
measurement uncertainty. The linear fits to the logarith- 
mic luminosities are plotted separately for the Seyfert 1 
galaxies, Seyfert 2 galaxies, and the combined sample, 
where we use the bisector of the fits of Lmir vs. Lx 
and Lx vs. Lmir- The coefficients of the fits are listed 
in Table [3] The linear fits to the Seyfert 1 and Seyfert 
2 galaxies are consistent with each other, indicating that 
the MIR emission is essentially independent of viewing 
angle and line of sight obscuration. The uncertainty of 
the fits to the limited Seyfert 1 sample alone is large, 
effectively a factor of 10 in the Lmir/Lx ratio. 

Considering the combined Seyfert 1 and 2 samples, 
both the diffraction-limited and PSF-fitting relationships 
suggest a dependence on luminosity, in the sense that 
Lx is reduced with greater MIR emission. We attribute 
this result to contamination in the MIR, with additional 
components of dust emission associated with star forma- 
tion present even on these small scales. Sources in ad- 
dition to the AGN do not affect the X-ray data we use, 
while they do contribute to the MIR flux. The galax- 
ies having excess MIR emission (the "contaminated" nu- 
clei) are preferentially those with higher MIR luminos- 
ity. Moreover, this effect is reduced when using the fixed 
physical aperture MIR measurement, in which the abso- 
lute amount of contamination does not vary significantly 
among the different galaxies. We can interpret this re- 
sult to mean that on scales of 100 pc, the star forma- 
tion in all these galaxies is comparable. The luminosity 
dependence in the diffraction-limited and PSF-fitting re- 
sults implies that the nuclear star formation (arising on 
even smaller scales) does vary. Star formation on such 
small scales in active galaxies has been directly mea- 
sured at near-IR wave lengths, showing a r ange of specific 
star formation rates (|Davies et al.1 12007) . Considering 
the MIR/X-ray relationships, while diffraction-limited or 
(especially) PSF-fitting measurements best reveal the re- 
processed AGN luminosity of an individual galaxy, sys- 
tematic differences introduce a luminosity dependence in 
the MIR/X-ray correlation that we do not attribute to 
properties of the AGN. 

We note that the sense of the result of I Gandhi et al.1 
(2009) is similar to ours, though weaker, with Lx ~ 
L°aiir° 06 ■ Thus, they do not claim any significant lumi- 




FlG. 1. — The MIR emission, which the reprocessing torus 
dominates, is strongly correlated with AGN bolometric luminosity, 
which X-ray luminosity (Lx) indicates. Here we show diffraction- 
limited MIR measurements. Solid blue, red, and black lines show 
the linear fits to the Seyfert 1 galaxies alone, the Seyfert 2 galax- 
ies alone, and all data combined, respectively. The uncertainty in 
the type 1 relationship is large, so the resulting fit is not signifi- 
cantly different from the type 2 or combined results. Observations 
of X-ray variability contribute to the effective error and are large 
for the Seyfert Is, while uncertainty of the intrinsic X-ray emission 
from the very absorbed AGNs accounts for the largest errors in the 
Seyfert 2s. Measurement uncertainty dominates the MIR errors. 
Additional data from the literature are plotted (blue diamonds); 
they do not contribute to the fits but are consistent with them. The 
filled symbol isNGC 4945 ( JTT}, and dotted lines show theoretical 
predictions ( jQO) . 




Fig. 2.— As in Figured] using PSF fitting of the MIR flux. The 
resulting fits are similar, and the correlation is again significant. 

nosity dependence in Lmir/Lx- While the luminosity 
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range of their sample is comparable, the average distance 
of their sources is greater. The luminosity effects we find 
are strongest in the measurements that best isolate the 
AGN torus, which may not be apparent at lower effective 

spatial resolut ion. 

In contrast. iMaiolino et al.l (|2007f ) found a luminosity 
dependence in the opposite sense, with MIR luminos- 
ity decreasing with increasing intrinsic AGN luminos- 
ity. That work investigated unobscured AGNs only, us- 
ing lower spatial resolution observations from the Spitzer 
Space Telescope. Concentrating on more distant and 
higher luminosity AGNs, including quasars, however, 
they attribute the dependence to genuine differences in 
the torus as a function of luminosity. Specifically, if 
high-luminosity quasars partially clear their immediate 
surroundings, having less dusty material available to re- 
process the intrinsic continuum would reduce the relative 
emergent MIR output. However, over the lower luminos- 
ity range we probe and utilizing high spatial resolution, 
contamination by star- heated dust is a stronger effect. 
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log 

Lmir [ er g s 1 (fixed physical aperture) 

Fig. 3. — As in FigurefT] where the MIR flux is measured within 
a fixed 100 pc aperture. The subsequent fits are again similar to 
the other results, and the correlation is again significant. 

If we extend the sample of Seyfert Is to include 
diffraction- limi ted measurements of nea rby galaxies from 
the literature (iHorst et al.ll2006l . I2008T I. the fitted func- 
tions remain consistent with each other over the differ- 
ent measurement techniques because of their large un- 
certainties. However, the strength and significance of the 
diffraction limited Seyfert 1 correlation is reduced, with 
p = 0.53 and P = 0.13. We plot these data in Figure [Q 
to show their overall agreement with our results, but the 
fitted relationships and statistical correlations of Table [3] 
do not use these measurements. 

3.2. Isotropic MIR Emission from Clumpy Models 

While a smooth optically and geometrically thick torus 
preserves the essence of AGN unification, accounting for 
viewing-angle dependent differences in the detectabil- 
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Wavelength (flm) 

Fig. 4. — Clumpy torus model spectra (Ncnkova et al.N2008bT) 
from different viewing angles show a range of silicate feature be- 
havior (around lO^tm) while remaining roughly isotropic through 
the MIR continuum. The inclination angles plotted are 0° (blue) 
through 30, 50, 60, 70, 80° to 90° (red). The spectra show very 
little variation from all type 1 (low inclination) views. The model 
parameters are a = 30° , No = 4, ry = 60, and Y = 30, with q = 1 
on the top, and q = 2 on the bottom. 

ity of spectrally broad optical emission lines, it can- 
not accommodate the observed MIR isotropy De- 
tailed radiative transfer models jP ier & Krolik 119921: 
i Granato fc Pan ose 1994; Efstat hTou &: Rowan-Robinson! 
Il995f ) all indicate significant differences in emergent MIR 
luminosity as a function of viewing angle for a given 
AGN luminosity. For any particular dust configuration, 
changes in viewing angle result in MIR luminosity dif- 
ferences of multiple orders of magnitude. Considering a 
range of simulated dust distributions, the resulting vari- 
ation is even larger. Fundamentally, the type 1 AGNs 
are preferentially brighter in the MIR because they offer 
unobscured views of the hot optically thin surface of the 
inner throat of the torus. In contrast, obscuration of the 
central engines of type 2 AGNs simultaneously hides this 
hottest (most emissive) dust behind a large optical depth 
of cool material. 

Clumpy torus models, in which the volume filling factor 
of the dusty material remains small, can readily account 
for the observed MIR isotropy while retaining the cen- 
tral features of AGN unification (e.g..lH5nig et alJl2006t 
iNenkova et~aH l2008ri ISchartmann et all l200gj l The 
principal difference is that these configurations can pro- 
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vide direct views of some hot and some cold cloud faces 
from all lines of sight, even when the active nucleus it- 
self is hidden. Thus, they reduce the sensitivity of MIR 
emission to viewing orientation, while simultaneously al- 
lowing for small tori and variability in line of sight ob- 
scuration. 

We u se the formalism and results of iNenkova et al.l 
( 2008a. b) to interpret these observations quantitatively. 
Individual clouds are optically thick in t he V band, given 
by tv , with the cool silicate dust of Osscnkop f et al.l 
(1992). The average number of clouds along an equa- 
torial ray is Nq, and they are distributed according to 
a power law in radius, oc R~ q , from the dust sublima- 
tion radius, Rd, to an outer radius R = YRd- The 
"torus" is concentrated in the equatorial plane and has 
a smooth edge, with a setting the distribution in eleva- 
tion; the average number of clouds along the line of sight 
Nlos = Aoexp(— (90° — i) 2 /a 2 ) at inclination angle i 
(measured from the pole). 

Figure [4] presents the resulting MIR spectral energy 
distributions as a function of inclination angle for some 
characteristic parameter values. They show extremely 
little variation over all type 1 views (i = to 30°). More 
importantly in the context of the present observations, 
the MIR continuum flux changes by only factors of a few 
from the extreme pole-on (i = 0°) to equatorial (i = 90°) 
lines of sight. One further consequence is that the sili- 
cate feature around 10/im can appear in emission (espe- 
cially from type 1 views) or absorption (especially from 
higher inclinations), but it is generall y weak, which is 
consistent with spectr al observations (|Hao et al.l 120071 : 
IThompson et al.ll2009f ). 

The emergent reprocessed MIR emission depends on 
the intrinsic AGN luminosity. We plot the theoreti- 
cal MIR/X-ray correlations as dotted lines in Figures Q] 
through [3] for the characteristic parameters of Figure |4] 
(adopting q = 2). These model relationships have no lu- 
minosity dependence (i.e., the slope of the lines is one), 
and show the offset factor of a few between the extreme 
pole-on (blue) and equatorial (red) viewing angles. The 
absolute scale depends on the bolometric correction to 
the X-ray luminosity. The plotted lines assume that the 
absorption-corrected 2-10 keV emission represents 5% of 
the intrinsic luminosity. The theoretical predictions are 
extremely close to the results of fitting the type 1 data 
alone for all MIR measurements. In contrast, most of 
the type 2 data lie to the right of the model lines, which 
further supports the interpretation that non-AGN MIR 
emission can be important even over very small angular 
scales. 

While the MIR emission is nearly isotropic, 
INenkova et al.l (|2008bl ) discuss several ways the clumpy 
model parameters affect the sensitivity of MIR emission 
to viewing angle. In general, the most compact tori 
produce the most isotropic emission. The compactness 
may be achieved by either reducing the outer torus 
extent (lowering the Y parameter), or by increasing the 
steepness of the radial power law, q, which concentrates 
the clouds toward the inner regions. Larger torus scale 
heights provide more isotropic clump distributions, 
so increasing a also reduces the sensitivity of MIR 
luminosity to viewing angle. Because increasing the 
number of clouds along radial rays (increasing Nq) 
approaches the smooth models, these distributions show 



greater MIR variation. The combination of No, cr, and i 
determines the number of clouds along the line of sight, 
with fewer clouds favoring greater variability and more 
views of directly-illuminated cloud surfaces. 

Here we emphasize the uniformity of the emergent MIR 
flux. We consider a range of parameter values, with 
N Q = 1-15, cr = 15-60, q = 0-3, and t v = 10-200 in 
over 11,000 models. The absolute 8.8/xm flux, Fg.g, does 
not vary by more than a factor of 600 considering all 
these models, whereas smooth torus calculations typi- 
cally show even larger variation within a single model re- 
alization. For any given clumpy configuration, i*g.g does 
not change by more than a factor of 35 as a function of 
viewing angle, and in more than 99% of the cases, the 
variation is less than a factor of 20. For typical values of 
the model parameters, changes in viewing angle do not 
result in MIR flux differences greater than a factor of 
5. Figure [5] illustrates these results, showing the flux at 
8.8/im as a function of viewing angle, scaled to the 8.8/im 
flux from the equatorial (edge-on) view of the cloud dis- 
tribution. Different panels show various values of a and 
q, and colors and line styles indicate Nq and Ty, respec- 
tively. As expected, isotropy at 8.8/im increases with 
increasing cr, increasing q, decreasing N q, and decreas- 
ing tv , consistent with the analysis of INenkova et al.l 
<|2008hV 

The efficiency of reprocessing intrinsic AGN emission 
to the MIR is sensitive to the total number of clouds 
in the torus, which iVo and cr govern. These parame- 
ters also determine the likelihood that the central en- 
gine is hidden. The less efficient reprocessors are more 
likely to be classified as type 1 AGNs, offering more un- 
obscured lines of sight. The effect of this selection bias 
is to reduce the net differences in MIR emission between 
AGN types. While for a given cloud distribution, type 
1 views show increased MIR flux over their type 2 coun- 
terparts, the distributions that are more likely to be ob- 
served as type 1 have reduced MIR flux relative to the 
cloud configurations that are more likely to be observed 
as type 2. Observational samples remain small, so con- 
clusions are not statistically robust, but they do suggest 
differences in the cloud distri butions between AGN types 
(jRamos Almeida et afll2009[ ). 

Observations indicate that typical values for parame- 
ters range over N ~ 5 -10, cr - 30-45°, and tv ~ 30-100 
(|Nenkova et al.l l2008bl and references therein). Specifi- 
cally c onsidering the 10 and 1 8 /mi silicate features to- 
gether, IThompson etaD (|2009h find N Q < 6, a = 45°, 
and Ty ~ 30-60. The MIR-emitting tori of the clumpy 
models are small. In terms of bolometric AGN lumi- 
nosity in uni ts o f 10 45 ergs _1 , L45, the outer radius 
R w 0.4Yy/L 45 pc, with Y — 30 or less. Even for 
the most luminous Seyfert galaxy of the sa mple, NGC 
3281 , R < 7 pc, assuming Lboi = 20Lx (j Elvis et al.1 
1994). The torus is not resolved in direct imaging, de- 
spite the good angular resolution of large ground-based 
telescopes. NGC 4151 offers the most favorable combina- 
tion of luminosity and proximity, yet here R corresponds 
to an angular scale < 0.07". This value is well below the 
diffraction limit of 0.27" at 8.8/im on an 8.1 m telescope. 
Intriguingly, resolving the torus may be possible on a 
30m telescope, where i?g.g = 0.074". MIR interferom- 
etry with MIDI at the VLTI currently offers sufficient 
resolution and has been used successfully to isolate AGN 
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tori tiJaffe et al.ll2004t iTristram et"aT1l2007h iRaban et all 
2009). However, this technique requires very bright tar- 
gets, so the number of AGNs available for study with 
this facility is li mited, with no mor e than a dozen being 
fe asible in tot a l (|Raban et alJl2008l). 

iHorst et all (|2008| ) and iGandhi et all (|2009f ) attempt 
to address possible contamination in their comparable 
measurements by considering their ability to resolve 
emission over larger scales. Specifically, they define 
sources to be "well-resolved" if the physical scale of an 
observation's angular resolution is less than 560Rd, which 
is many times larger than the outer torus radius R . 
The clumpy model predictions and basic considerations 
of the temperature of dust that emits significantly in the 
MIR argue instead that the relevant size scales are much 
smaller, so "resolution" on these scales of 560Rd does not 
limit the measurements to the torus alone. However, the 
empirical support these authors find for the utility of this 
resolution limit may be a consequence of dust heated by 
stars and star formation, which can dominate on scales 
larger than their 100 pc resolution limit (w 560i?d for 
their typical Lx = 10 43 ergs _1 ). We find these non- 
AGN contributions to be significant on all scales. The 
flux in the fixed physical aperture is (on median) 20% 
higher than the PSF-fitted value. On scales smaller than 
100 pc, the nuclear star formation produces the observed 
luminosity dependence in the MIR/X-ray relationship. 

One completely different explanation for the observed 
isotropy is that the t orus i s simply optically thin in the 
MIR, as lHorst et al.l (|2006| ) considered. However, an op- 
tically thin torus would show 10/im silicate emission in 
all cases, contrary to spectral observations. While silicate 
absorption attributable to the torus is generally weak, it 
is typical of type 2 AGNs, which requires optically thick 
material along the line of sight. (A smaller column in 
a cold foreground screen could alternatively account for 
the absorption, but this could not simultaneously be the 
source of the IR emission.) Furthermore, the Compton 
thick AGNs require total ry > 800 along the line of sight 
(assuming standard gas-to-dust ratios) . This value corre- 
sponds to total r > 14 at 8.8um, although dust-free gas 
(located within the dust sublimation radius) may also 
contribute to their X-ray absorption. 

4. CONCLUSIONS 

Intrinsic X-ray and reprocessed MIR luminosity of the 
active nuclei of Seyfert galaxies are strongly correlated. 
Emphasizing measurements at high spatial resolution, we 
find a luminosity dependence on the correlation over the 
range we probe when using two MIR measurement tech- 
niques that help to isolate the central unresolved emis- 
sion. We attribute this effect to the presence of addi- 
tional sources in the MIR, such as by emission by dust 
that young stars heat. Because this result emerges only 
using the two smallest scale measurements, it implies 
that significant nuclear star formation is present even 
on very small (sub-100 pc) spatial scales in these AGNs. 
However, we find no luminosity dependence in the MIR 
fluxes measured in fixed 100 pc apertures, suggesting 
that the stellar contributions on these scales are com- 
parable in all galaxies. 

We find no significant difference in the correlations be- 
tween Seyfert 1 and Seyfert 2 galaxies, given the large (~ 
factor of 10) uncertainty in the Seyfert 1 Lmir/ Lx ratio. 



Thus, we conclude that the MIR emission of the AGN 
torus is nearly isotropic on all the small scales we probe 
here. We interpret these results as support for clumpy 
torus models, in which the MIR emission is insensitive to 
viewing angle because both illuminated and dark cloud 
faces contribute in part to observed flux over all viewing 
angles. In contrast, smooth obscuration models produce 
significantly stronger (by several orders of magnitude) 
MIR emission from type 1 AGNs, which exclusively al- 
low direct views of the hottest dust surfaces. 

We use numerical radiative transfer calculations to 
demonstrate explicitly the small differences in the spec- 
tral energy distributions for extreme pole-on (type 1) and 
edge-on (type 2) views of a clumpy torus. For typical 
model parameters, the MIR flux changes by less than a 
factor of 5 over all viewing angles. Considering a wide 
range of different cloud distributions, the observed varia- 
tion is nearly always less than a factor of 20. Several fac- 
tors that contribute to increasing the resulting isotropy 
include increasing the spherical symmetry of the cloud 
distribution, reducing the number of clouds along radial 
rays, and making the clumpy torus more compact. Fur- 
thermore, flux differences as a function of viewing angle 
decline at longer wavelengths. 

While the clumpy models predict nearly isotropic MIR 
emission, selection effects may further contribute to min- 
imizing the observed difference between AGN types, 
whose classification is sensitive to viewing angle. The tori 
of the AGNs that are most likely to be observed with- 
out obscuration reprocess the intrinsic emission less effi- 
ciently than those that contain more dusty clouds (which 
are more likely to be classified as type 2s). Consequently, 
the relative MIR flux of observed type Is may be prefer- 
entially suppressed, while the observed type 2s may be 
preferentially MIR-enhanced. 

An optically thin torus could account for the MIR 
isotropy, but it could not produce the observed spec- 
tral differences. Specifically, silicate absorption is typical 
of Seyfert 2s, and weak silicate features (in emission or 
absorption) are characteristic of Seyfert Is. The clumpy 
torus is not optically thin in the MIR, which allows for 
these spectroscopic differences. Because both directly- 
illuminated and dark cloud faces contribute to the net 
observed spectrum, the resulting silicate features are al- 
ways weak. 
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Fig. 5. — MIR flux variation as a function of viewing angle. Each curve covers one clumpy dust distribution, and the flux is scaled to 
the value observed from the equatorial (edge-on) view. These results illustrate that the emergent 8.8/im flux does not vary by more than 
factors of about 5 over all viewing angles, for typical parameter values (Nq < 6, Ty < 60). Making the torus effectively more compact 
(increasing q) or making the cloud distribution more spherically symmetric (increasing <r) increase the MIR isotropy, as do decreasing No 
and tv- 



TABLE 1 
Observation Log 



Galaxy 


Instrument 


Telescope 


Observation Date 


Filter 


On-Source Time (s) a 


Circinus 


T-ReCS 


Gemini S 


2004 Feb 01 


Si2 


565 


IC 5063 


T-ReCS 


Gemini S 


2005 Jul 18 


Si2 


130 


NGC 1068 b 


Michelle 


Gemini N 


2004 Aug 10 


Si5 


4 








2005 Jan 31 


Si5 


4 


NGC 1365 


OSCIR 


CTIO 


1998 Feb 12 


N 


482 


NGC 1386 


T-ReCS 


Gemini S 


2003 Dec 06 


N 


217 


NGC 1566 


T-ReCS 


Gemini S 


2005 Sep 18 


Si2 


152 


NGC 2992 


Michelle 


Gemini N 


2006 May 05 


N' 


730 


NGC 3081 


T-ReCS 


Gemini S 


2006 Jan 25 


Si2 


130 


NGC 3227 


Michelle 


Gemini N 


2006 Apr 07 


N' 


300 


NGC 3281 


T-ReCS 


Gemini S 


2004 Jan 30 


N 


261 


NGC 4151 


OSCIR 


Gemini N 


2001 May 07 


N 


360 


NGC 4388 


Michelle 


Gemini N 


2006 May 12 


N' 


549 


NGC 4945 


T-ReCS 


Gemini S 


2006 Mar 17 


Si2 


261 


NGC 5506 


Michelle 


Gemini N 


2006 Apr 06 


N' 


546 


NGC 5728 


T-ReCS 


Gemini S 


2005 Jul 08 


Si2 


130 


NGC 7172 


T-ReCS 


Gemini S 


2004 May 13 


N 


305 


NGC 7582 


OSCIR 


Gemini S 


2001 Dec 13 


N 


608 



a In guided beam only. b Both observations were used for flux calibration, with spatial profile 
measurements of only the 2004 data, obtained under better seeing conditions. 
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TABLE 2 
Nearby Seyfert Galaxy Sample 



Galaxy 


z 


D 


Type L x a 


Tdl b 
MIR 


fF/ c 
MIR. 


jphys d 
MIR 


PSF Scale X-ray Reference(s) 






(Mpc) 


log (erg s" 1 ) 


log (erg s- 1 ) 


log (erg s 


log (erg s 1 ) 


(%) 



42.61 


42.57 


42.75 


90 


1 


43.57 


43.57 


43.61 


100 


2 


44.12 


43.81 


44.34 


50 


3, 4 


42.82 


42.77 


42.83 


90 


5 


42.04 


41.88 


42.22 


70 


3 


41.73 


41.73 


41.77 


100 


6 


42.82 


42.78 


42.85 


90 


7 


42.58 


42.58 


42.60 


100 


6 


42.53 


42.53 


42.50 


100 


8, 9 


43.49 


43.40 


43.61 


80 


10 


43.00 


42.95 


43.08 


90 


11, 12, 13 


42.94 


42.90 


42.95 


90 


14, 15 


40.25 


39.25 


41.83 


10 


16, 17 


43.28 


43.28 


43.28 


100 


18, 9 


42.20 


42.20 


42.20 


100 


6 


42.55 


42.51 


42.58 


90 


19, 20, 21 


42.32 


42.23 


42.49 


80 


22 



Circinus 
IC 5063 
NGC 1068 
NGC 1365 
NGC 1386 
NGC 1566 
NGC 2992 
NGC 3081 
NGC 3227 
NGC 3281 
NGC 4151 
NGC 4388 
NGC 4945 
NGC 5506 
NGC 5728 
NGC 7172 
NGC 7582 



0.0014 
0.0113 
0.0038 
0.0055 
0.0029 
0.0050 
0.0077 
0.0080 
0.0039 
0.0107 
0.0033 
0.0084 
0.0019 
0.0062 
0.0094 
0.0087 
0.0053 



Additional galaxies' 3 
MCG -6-30-15 0.0077 
NGC 3783 0.0097 
NGC 4593 0.0090 
NGC 7314 0.0050 



4.0 
47.9 
14.4 
23.0 
12.2 
21.2 
32.6 
33.7 
16.3 
45.1 
14.0 
35.5 

3.7 
26.1 
39.5 
36.7 
22.2 

32.7 
41.1 
38.0 
21.1 



2.0 
2.0 
2.0 
1.8 
2.0 
1.0 
2.0 
2.0 
1.5 
2.0 
1.5 
2.0 
2.0 
1.9 
2.0 
2.0 
2.0 

1.2 
1.5 
1.0 
1.9 



41 76+ 008 
—0.10 

42.93t°;°| 
42.78jll;gl 
42.18l°;« 



41.81 



+0.08 
0.09 
r n +0.13 



41.501 . lg 
42.181°;^ 

42 n +0.08 
' 1 -0.10 

43.18+ - 08 



-0.09 
-0.04 
-0.79 

42.881°-^ 
41.331°^ 
43.06l°;« 

,+0.07 
-0.10 
-+0.34 
-0.24 

42.691°;?? 



42.94 



42.21] 
42.75" 1 



42.59 



+ 0.21 
0.15 



43.15t° £ 
42.88lO.og 

-,+0.30 
-0.29 



42.22 i 



43.09 
42.55 
42.41 
41.69 



23 
24, 25 
26 
27 



References. — (1) Soldi ct al. 2005; (2) Turner ct al. 1997; (3) Lcvcnson ct al. 2006; (4) Iwasawa et al. 1997; (5) Risaliti et al. 2005; (6) This 
work: (7) Yaqoob ct al. 2007; (8) Lamer ct al. 2003; (9) Uttlcy & McHardy 2005; (10) Vignali fc Comastri 2002; (11) Bcckmann ct al. 2005; (12) 
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al. 2000; (18) Lamer et al. 2000; (19) Awaki ct al. 2006; (20) Akylas ct al. 2001; (21) Smith & Done 1996; (22) Turner ct al. 2000; (23) Miniutti 
ct al. 2007; (24) Do Rosa et al. 2002; (25) Kaspi ct al. 2002; (26) Reynolds et al. 2004; (27) Yaqoob ct al. 2006. 

a Absorption-corrected 2-10 kcV X-ray luminosity. 15 Diffraction-limited aperture measurements Measurement from PSF fitting. d Fixed physical 
(100 pc) aperture measurement. 6 From Horst et al. (2008). 



TABLE 3 
Correlation Results 



Sub-Sample 




Diffraction Limited 






PSF Fitting 






Physical Aperture 




P 


P a 


b 


P 


P a 


b 


P 


P a 


b 


All galaxies 
Seyfert 1 
Seyfert 2 


0.78 
0.90 
0.74 


1.8c-03 0.76 ± 0.06 
7.2e-02 1.0 ±0.2 
1.4e-02 0.70 ± 0.07 


9.9 
-2.0 
12.6 


0.79 
0.90 
0.78 


1.5e-03 0.68 ±0.06 
7.2c-02 1.1 ±0.2 
l.Oc-02 0.61 ±0.06 


13.3 
-2.7 
16.4 


0.77 
0.90 
0.73 


2.0c-03 0.93 ± 0.06 
7.2e-02 1.0 ±0.2 
1.6e-02 0.92 ±0.07 


2.7 
-2.1 
3.0 



Note. — Spearman rank correlation coefficient, p, and two-sided probability, P. Linear fit parameters, for \og(Lx) — alog(Z/MJR) + 



